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I.  INTRODUCTION 


A.  Background 

A  theoretical  investigation  of  the  phase  and  amplitude  behavior  of  light 
pulses  amplified  in  a  GaAs  double  heterostructure  device  including  the  effect 
of  intraband  relaxation  processes  applied  to  tbe  individual  momentum  states  is 
presented.  Fast  intraband  processes  were  considered  in  early  work  concerning 
the  quantum  mechanical  treatment  of  fluctuations,  relaxation,  and  noise  in 
semiconductor  lasers  [1,2].  These  processes  were  assumed  to  be  so  fast  com¬ 
pared  to  all  times  of  interest  that  each  electron  and  corresponding  hole  was 
considered  to  be  in  thermal  equilibrium  with  all  the  other  electrons  and 
holes.  In  this  report  we  do  not  require  that  processes  be  fast,  and  we  calcu¬ 
late  the  transition  for  each  state  separately.  These  processes  were  also 
considered  in  a  nonlinear  treatment  for  laser  diodes  [3],  modeled  in  terms  of 
relaxation  times  required  for  electrons  and  holes  as  an  entirety  to  come  to 
quasi  equilibrium.  The  concept  of  very  fast  intraband  relaxation  processes 
resulted  in  the  use  of  a  single  two-level  system  representing  all  the  elec¬ 
trons  and  holes.  With  this  concept  intraband  relaxation  processes  could  be 
Ignored  and  a  single  rate  equation  written  for  the  population  inversion.  This 
equation,  plus  a  rate  equation  for  photon  number  density  and  the  Fabry-Perot 
(F-P)  oscillator  condition  for  simple  resonators,  permitted  analysis  of  the 
semiconductor  laser  applied  to  communication  problems  [4,5].  An  extensive 
theoretical  analysis  of  the  gain,  frequency  bandwidth,  and  saturation  output 


power  o f  double  heterostructure  laser  amplifiers  was  made  possible  because  of 
the  inclusion  ot  multimode  rate  equations  16 J.  More  recently,  an  interest  in 
hign-data-rate  systems  and  a  requirement  to  determine  the  cross  interference 
of  various  frequency  modes  in  the  lasers  has  stimulated  a  reconsideration  of 
intraband  relaxation  processes.  fi.  Yamada  and  Y.  Suemetsu  [7]  considered 
these  processes  in  a  density  matrix  formulation  of  the  problem  of  gain  sup¬ 
pression  in  undoped  injection  lasers.  A  similar  analysis  was  performed  by 
K.  Kazarinov,  C.  Henry,  and  R.  Logan  [ 8 J . 

With  the  advent  of  techniques  used  to  suppress  regeneration,  interest  has 
also  risen  in  the  traveling  wave  amplifier  (TWa).  Luk'yanov  et  al.  [9]  ob¬ 
tained  experimental  results  on  a  TWA  by  orienting  a  stripe  homojunction  laser 
at  an  angle  to  the  cleaved  end  faces  of  the  crystal.  An  analysis  was  made  of 
the  amplifier  as  a  function  of  the  pumping  rate  and  input  signal.  Small  sig- 
nal  gains  exceeding  10  were  reported.  Otsuka  [10]  proposed  and  theoretically 
investigated  laser-amplifier-based  integrated  optical  circuits  by  means  of 
simplified  TWA  rate  equations  without  intraband  processes  included. 

In  this  report  the  theory  of  lesevgi  and  Lamb  [11]  is  applied  to  the 
semiconductor  TWA;  it  differs  from  previous  work  by  retaining  the  effect  of 
the  Intraband  relaxation  on  the  gain.  In  this  report,  for  the  first  time, 
phase  has  been  considered  in  an  analysis  of  hetero junction  amplifiers.  The 
results  of  this  study  are  important  to  the  design  of  amplifiers  that  must 
produce  an  output  in  which  phase  is  fixed  with  respect  to  the  input  phase  and 
few  multiple-pulse  effects  are  exhibited.  These  amplifiers  may  be  useful  in 
heterodyne  or  homodyne  communication  systems,  as  components  of  a  phase  locked 
oscillator 'amplifl'  high  power  laser  source,  and  for  optical  signal 


processing 
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B.  Summary  of  Results 

The  phase  of  the  output  relative  to  the  input  is  strongly  dependent  on 
frequency,  the  intraband  collision  induced  homogeneous  broadening,  and  the 
instantaneous  density  of  states  of  the  semiconductor.  It  is  indeed  possible 
to  obtain  a  zero  phase  change  for  the  amplified  signal.  Small  signal  gains 
greater  than  leP  are  predicted,  in  agreement  with  Ref.  9.  Saturation  of 
intensity  will  occur  as  a  consequence  of  the  model.  The  effect  of  current 
limiting.  Introduced  by  a  resistance  (occurring  naturally  or  introduced  to 
prevent  diode  burnout)  in  series  with  the  junction,  is  demonstrated  to  reduce 
the  gain  of  the  amplifier.  For  pulsed  operation,  the  amplifier  introduces 
pulse  compression,  as  expected,  and  also  Introduces  a  nonuniform  phase  varia¬ 
tion  with  time.  This  variation  may  limit  phase  coherency  of  the  output  with 
respect  to  the  input  field.  Multi-pulse  effects  will  be  exhibited  for  pulse 
separations  shorter  than  several  pumping  times. 

C.  Theoretical  Assumptions 

Following  Icsevgi  and  Lamb  [  10 ]  a  self-consistent  semiclasslcal  descrip¬ 
tion  of  the  interaction  between  matter  and  radiation  field  is  used  in  which 
the  medium  is  treated  quantum  mechanically  and  the  radiation  field  is  treated 
according  to  Maxwell's  theory.  The  investigation  is  carried  out  in  this 
manner  because,  in  contrast  to  the  approaches  of  Refs.  4,  5,  6,  8  and  9,  the 
phase  can  be  included.  Saturation  occurs  as  a  consequence  of  the  model,  and 
the  intraband  relaxation  processes  can  be  introduced  in  a  more  physical  way. 
The  starting  point  of  the  analysis  is  the  density  matrix  equations,  which 
Include  the  phenomenological  rate  constants.  The  sources  of  various  terms  in 
the  equations  are  discussed. 


Self-consistency  is  invoked  by  requiring  that  the  electric  field  propa¬ 
gating  in  the  medium  introduce  a  macroscopic  polarization  which,  in  turn,  acts 
as  a  source  term  to  Maxwell's  equations.  The  following  assumptions  are  made 
to  simplify  the  problem. 

1.  The  electromagnetic  field  E(z,t)  is  a  scalar  representing  a  uniform  plane 
wave  polarized  in  the  x  direction  and  propagating  in  the  z  direction. 
Questions  concerning  regeneration  and  modal  structure,  therefore,  are 
ignored  and  a  TWA  analysis  is  performed. 

2.  Amplification  in  a  standard  double  heterostructure  device  is  con¬ 
sidered.  The  medium  consists  of  a  group  of  two-level  electron  states 
generated  by  applying  a  bias  voltage  across  a  P-n-N  junction.  The  con¬ 
duction  and  valence  band  occupations  are  described  by  the  Fermi-Dirac 
distribution  functions.  Electron  transitions  between  two  energy  level 
distributions,  rather  than  between  two  discrete  levels,  must  be  con¬ 
sidered  for  interband  absorption  or  emission  of  radiation  at  a  given 
frequency.  To  further  simplify  the  problem,  the  parabolic  band  approx¬ 
imation  is  made  and  it  is  assumed  that  the  k-selectlon  rule  is  valid. - 
This  choice  will  not  significantly  affect  the  spontaneous  emission 
spectrum  for  energies  greater  than  the  band  gap  energy,  which  establishes 
the  relationship  between  the  circular  frequency  u  of  a  given  transition 
and  its  corresponding  wavenumber  k  as 


*2.2 

"Kw  »  E  +■  — x — 
g  2m 


where  Eg  is  the  band-gap  energy,  mr  is  the  reduced  mass,  and  "h  is 
Planck's  constant  divided  by  2ir.  Amplification  will  generally  occur  in  a 


frequency  band  bounded  by  Che  band-gap  and  the  difference  between  quasi- 
Fermi  levels  of  the  conduction  (Fc)  and  valence  (Fy)  bands. 

E  <  "Ku>  <  F  -  F  (2) 

g  c  v 

It  will  be  seen  later,  for  steady  state,  that  the  gain  and  phase  of  the 
wave  is  a  convolution  integral  of  a  collision  line  profile  and  a  density 
of  single  particle  states.  For  this  reason  the  parabolic  and  k-selection 
rule  assumptions  are  not  as  restrictive  as  it  first  seems  since  other 
models  for  the  ,£pnd  structure,  such  as  those  postulated  by  Laser  and 
Stern  111],  could  easily  be  substituted  into  the  final  results.  On  the 
other  hand,  the  model  postulated  by  Aleksanian,  et  al.  [ 1 2 J  could  not  be 
easily  used,  because  single  particle  states  are  not  assumed  for  that 
analysis. 

The  dipole  approximation  holds  for  the  interaction  between  the  radiation 
field  and  the  electronic  states. 

The  two-level  states  are  coupled  only  through  their  interaction  with  the 
common  electromagnetic  field  E(z,t). 

Phenomenological  rate  constants  are  included  to  describe  the  effects  of 
pumping,  spontaneous  emission,  and  phonon-electron  or  electron-electron 
collision  Induced  intraband  decay.  Interband  collision-induced  decay 
is  not  included,  because  it  would  make  the  problem  only  slightly  more 
difficult  and  not  reveal  significantly  more  information. 

The  optical  frequency  u  representing  a  given  transition  is  larger  than 
the  optical  linewidth. 
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The  rotating  wave  approximation  applies,  i.e.,  terms  involving  the  crea¬ 
tion  of  a  photon  and  simultaneous  creation  of  the  upper  state  leading  to 
harmonics  of  the  optical  field  are  ignored.  Optical  harmonic  generation 
is,  therefore,  excluded. 

8.  The  slowly  varying  wave  approximation  is  invoked. 

9.  All  pulses  of  interest  are  considerably  longer  than  the  intraband  relaxa¬ 
tion  time. 

In  Section  II,  equations  describing  the  phase  and  amplitude  behavior  of 
waves  propagating  in  the  heterojunction  device  are  given.  Numerical  results 
are  presented  in  Section  III,  and  a  brief  summary  appears  in  Section  IV. 

The  appendix  contains  a  more  detailed  derivation  of  the  theory  presented 
in  the  main  body  of  the  text.  Note  that  there  is  a  notational  difference  from 
the  main  text.  In  particular,  subscripts  2  and  1  of  the  text  are  c  and  v  in 
the  appendix.  The  discussion  pertaining  to  the  calculation  of  a  finite  phase 
change  on  amplification  is  not  included  in  the  appendix. 


II.  THEORY 


A.  General  Equations 

The  density  matrix  is  introduced  to  represent  the  state  of  the  two  levels 
assumed  for  the  medium.  Following  Icsevgi  and  Lamb  [10],  the  equations  of 
motion  in  component  form  for  the  density  matrix  are  derived  using 
Schrodinger 's  equations  and  the  dipole  approximation  as 


3p. 
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(3) 


(4; 

(5) 


P21 


(6) 


where  P2  is  the  density  of  upper  states  at  time  t,  macroscopic  position 
z  [6] ,  and  energy  level  or  momentum  k.  The  upper  level  represents  a  minority 
carrier  electron  in  the  conduction  band  and  a  corresponding  hole  in  the 
valence  band.  The  density  of  lower  states,  corresponding  to  an  empty  state  in 
the  conduction  band  and  an  electron  in  the  valence  band,  is  represented  by  p^. 
The  electromagnetic  field  is  represented  by  E.  The  contribution  to  the  total 
polarization  caused  by  a  given  k  state  produced  by  the  above  density  of  states 
is  given  by  the  partial  polarization 
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P(k, z,t )  *  p[p21(k,z,t)  +  c.c.] 

where  p  is  the  expectation  value  of  the  dipole  moment. 

The  other  terms  in  Eqs.  (3)  through  (6)y  ,Y,Y*P«*  and  p,  have 

zi  s  p  zo  lo 

been  introduced  phenomenologically  as  indicated  in  assumption  5.  The 
dephasing  rate  Y21  is  the  inverse  of  the  dipole  moment  relaxation  time  T2« 
This  rate  is  determined  by  the  intraband  relaxation  processes  including  polar- 
optical  -phonon  inelastic  scattering)  ionized  impurity  and  heavy-hole  elastic 
scattering,  and  electron-electron  inelastic  scattering  [13].  At  high  doping 
levels  and  high  levels  of  injected  carriers  the  dominant  scattering  process  is 
electron-electron.  Because  this  process  is  relatively  insensitive  to  energy 
level  [14],  the  dephasing  rate  is  assumed  independent  of  k  in  this  report.  We 
note  that  the  introduction  of  the  intraband  relaxation  rate  by  means  of  Eqs. 
(3)  through  (5)  is  equivalent  to  the  assumption  of  an  exponential  decay  to 
equilibrium.  In  Section  II. B,  this  assumption  is  modified  to  obtain  a  finite 
phase. 

Steady-state  pumping  of  the  amplifier  is  accomplished  by  maintaining  « 
constant  bias  across  the  semiconductor  junction.  Individual  k  state  pumping 
is  a  complicated  process  [1]  Involving  intraband  relaxation  processes,  dif¬ 
fusion  times  for  the  carriers,  and  the  instantaneous  voltage  across  the  junc¬ 
tion.  During  amplification  the  population  levels  are  reduced  by  stimulated 
emission.  This  causes  a  reduction  in  the  instantaneous  voltage  across  the 
junction.  Because  this  voltage  determines  the  quasi-Fermi  levels,  the  fast 
intraband  relaxation  processes  tend  to  restore  the  population  to  the  instan¬ 
taneous  values  determined  by  these  Fermi  levels.  The  Intrinsic  diode  capaci- 


16 


tance,  resistance,  and  carrier  mobility  as  well  as  external  circuit  parameters 

determine  the  way  in  which  these  levels  return  to  the  equilibrium  values  in 

the  absence  of  an  electromagnetic  optical  pulse.  The  first  attempts  to 

accurately  describe  the  electrical  characteristics  of  the  behavior  of  a  diode 

have  only  recently  been  published  [15,16].  A  mathematical  description  of  the 

processes  would  result  in  an  additional  set  of  equations  to  be  solved 

simultaneously  with  Eqs.  (3)  through  (6).  However,  because  the  optical  pulse 

behavior  rather  than  the  detailed  pumping  behavior  is  of  Interest,  a  single 

pumping  rate  constant  determined  by  the  slowest  process  is  used  117]. 

The  pumping  process  discussed  above  restores  the  density  of  states  to 

equilibrium  values,  p«  and  p.  ,  determined  by  the  Fermi  levels  in  the  absence 

zo  lo 

of  stimulated  or  spontaneous  emission.  These  values  are 


fcC1  *  V 
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where  fc  and  fv  are  the  Fermi  distribution  functions  given  by 


The  quantities  mc  and  mv  are  the  electron  and  hole  masses,  respectively,  Fc 

and  Fv  are  the  quasi-Fermi  levels,  K  is  Boltzmann's  constant,  and  T  the 

absolute  temperature.  The  final  constant  expressed  in  Eqs.  (3)  and  (4)  is 

the  spontaneous  emission  rate  Y  . 

s 

The  electromagnetic  field  equation  corresponding  to  the  above  set  of 
equations  determining  the  space-time  behavior  of  the  two-level  states  is 


+  p  a 
o 


2  -2  a‘E 

+  n  c  — IT 

°  3t 


P(k,z,t) 


where  the  conductivity  o  accounts  for  any  linear  losses  incurred  in  the 

absorbing  medium,  c  is  the  velocity  of  light  in  vacuum,  and  p  is  the  magnetic 

o 

permeability.  Note  that  the  total  polarization  P(z,t)  enters  Eq.  (12)  as  the 
sum  of  partial  polarizations  P(k,z,t),  whereas  Eqs.  (1)  through  (4)  can  be 
written  in  terms  of  the  partial  polarization  of  a  single  k  state.  The  intrin¬ 
sic  polarization  of  the  material  resulting  from  the  background  crystal  matrix 

is  included  in  the  index  of  refraction  n^. 

o 

Simplification  of  Eqs.  (3)  through  (6)  and  (12)  results  from  assumptions 
6  through  8  in  Section  I  and  the  use  of  retarded  time.  Assumption  8  means 
that  the  electromagnetic  field  can  be  written  in  terms  of  a  slowly  varying 
amplitude  and  phase  4  as 


E(z,t)  ■  <£(z,t)  cos[ vt  -  nz  +t(z,t)] 


p..-  .... 
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while  the  corresponding  partial  polarization  may  be  written  in  terms  ol 
in-phase  and  out-of-phase  components  C  and  S,  respectively,  as 


P(k,z,t)  ■  C(k,z,t)  cos(vt  -  nz  +  $)  +  S(k,z,t)  sin(vt  -  nz  +  $)  (14) 


where  higher  harmonics  of  the  field  are  neglected.  The  circular  frequency  v 
has  a  nominal  value  close  to  the  center  of  the  bandwidth  of  the  input  pulse, 
whereas  n  is  the  wavenumber  given  by 


Finally,  by  using  retarded  time  t,  2',  and  assumption  9,  a  coupled  set  of 
simultaneous  equations  can  be  written  for  the  electromagnetic  field  amplitude 
$  and  phase  4>  as  well  as  for  the  upper  and  lower  level  density  of  states  p^ 
and  p^. 


||  -  -L<f  +  2g k2(p2  -  Pj)  <^(u.)dk 


9  U  - 


— #C  k2(u>  -  v)(p2  -  Pj)  if(a))dk 


rm  "Yp(p2  '  P2o)  ‘  V2  +  ft 


J7~  “  'Yp(pl  “  Plo)  +  Ysp2~  7K 


The  Inclusion  of  2ir  in  v  is  from  Ref.  10. 
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where  e  Is  Che  permitivity  of  free  space  and  34>/3t  has  been  neglected  compared 
to  u>  -  v.  Note  that  approximation  9  means  that  the  polarization  is  “locked" 
to  the  field.  As  soon  as  the  pulse  has  passed  a  given  location,  the  polar¬ 
ization  at  that  location  also  drops  to  zero  and  the  remaining  inversion  will 
not  add  energy  into  the  tail  of  the  pulse.  This  means  that  phenomena  such 
aa  pulse  ringing,  as  discussed  by  lesevgi  and  Lamb  [10]  and  Allen  and 
Eberly  [18],  are  excluded.  Specialization  to  steady-state  and  pulsed  regimes 
is  considered  next. 

B.  Steady  State 

Steady-state  operation  of  the  amplifier  implies  setting  all  time  deriva¬ 
tives  in  Eqs.  (16)  through  (23)  equal  to  zero.  This  set  of  equations  then 
reduces  to 


-  -2L  I(z’)  +  2gnI(z’)  / dk  k2  #(u>)  N(k,z ’ ) 


~  f  dk  k2(u>  -  v)  ^Cw)  N(k.z’) 
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where  the  dimensionless  intensity  I(z')  has  been  introduced  and  is  defined 
here  as 


I  - 


The  population  difference  N(k,z*) 


+  V 


Pj  is  given  by 


( 26 ) 


N(k,2') 


l(yl>-y,)/<yJ>*y,)l  p2o  -  p1o 

i  +  &(  <0)1 


(27) 


Note  that  saturation  effects  are  given  implicitly  because  N,  appearing  in  the 
effective  gain  in  Eq.  (24),  has  the  intensity  as  an  inverse  factor.  Other 
steady-state  theories  include  the  effect  of  saturation;  however,  this  effect 
was  previously  introduced  phenomenologically  and  did  not  exhibit  an  integral 
dependence  on  momentum,  which  occurs  here  as  a  result  of  presence  of  the 
intraband  relaxation  processes*  The  equation  for  phase  for  the  he tero junction 
amplifier  fEq.  (23)]  is  original  in  this  report. 

The  phase  of  the  electromagnetic  field  is  an  expression  of  the  real  part 
of  the  index  of  refraction  and  as  such  is  known  to  be  finite.  On  the  other 
hand,  the  phase  as  expressed  by  Eq.  (25)  is  unbounded  with  the  use  of  the 
parabolic  approximation  and  an  infinite  upper  limit  on  k.  Although  most 
authors,  in  discussing  the  dispersion  relations  for  solids,  give  an  Integral 
expression  for  the  real  part  of  the  Index  of  refraction,  they  do  not  evaluate 
it  11-4-21],  but  use  experimentally  determined  values  Instead.  We  evaluate  the 
integral  appearing  In  Eq.  (25)  to  determine  the  phase  both  by  multiplying  the 
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Lorentzian  profile  lEq.  (21)]  by  a  Gaussian  profile  to  cut  off  the  quadrat- 
ically  decaying  tails  of  the  Lorentzian  and  by  setting  an  upper  limit  on 
momentum  (k)  integration.  A  previous  inclusion  of  intraband  relaxation 
processes  in  the  calculation  of  gain  [22]  was  based  on  many  body  theory  [12J 
and  the  spectral  weight  function  [24].  The  Lorentzian  spectral  weight 
function  used  in  Ref.  23  was  derived  assuming  an  exponential  decay,  which  is 
consistent  with  Eqs.  (3)  through  (5).  This  approximation  is  perfectly 
reasonable  for  a  gain  calculation  but  results  in  unbounded  phase.  Justifica¬ 
tion  for  the  modification  to  the  Lorentzian  profile  is  based  on  the  statisti¬ 
cal  theory  of  pressure  broadening  [24J  and  the  use  of  a  coulomb-screened 
potential  [23]  for  electron-electron  scattering.  The  use  of  the  screened 
potential  in  the  statistical  theory  of  line  broadening  results  in  Gaussian 
decay  in  the  wings  of  the  line.  The  line  center,  however,  decays  slightly 
less  rapidly  tnan  the  Lorentzian,  and  the  switcn  to  Gaussian  decay  occurs 
between  10  and  40  line  widths.  The  spectral  line  profile  that  was  used  in  the 
calculation  was  then  chosen  as  a  Lorentzian  multiplied  by  a  Gaussian  with  a 
line  width  (at  e-1)  of  20 

Examination  of  Bloom  and  hargenau’s  quantum  theory  of  spectral  line 
broadening  [26]  reveals  that  the  optical  field  causes  transitions  between 
"collision-smeared"  states.  The  amount  of  smearing  is  obviously  energy 
dependent.  If  there  are  no  free  electrons  above  a  given  energy  level,  smear¬ 
ing  beyond  that  energy  is  not  possible.  The  Fermi  distribution  for  electrons 
in  the  conduction  band  effectively  determines  the  maximum  energy  electrons 
that  can  occur  and  thereby  determines  the  upper  level  on  the  k-integration. 
This  upper  limit  was  chosen  at  9.1  x  10^  cm  *,  which  yields  a  value  of  0.01 


22 


for  the  conduction-band  Fermi  function.  It  is  interesting  to  note  tnat  this 
energy  cutoff  is  consistent  with  the  potential  energy  of  electrons  separated 

1  O 

by  a  screening  distance  of  30  A  given  for  an  electron  density  of  10  1271. 
This  was  the  screening  distance  used  in  determining  the  Gaussian  profile 
discussed  above. 

To  permit  numerical  predictions,  the  dipole  moment  p  must  be  determined. 
From  Ref.  28  we  can  write 
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2n  Kc  n  c  iy 

r  °  o  ]  <2 
L  t/2  J 
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(28) 


where  K  is  an  absorption  parameter  experimentally  measured  as  6000  cm  1 
(eV)-^/^.  For  values  of  reduced  mass  and  the  transition  energy  given  in 
Ref.  28  and  with  nQ  *  3.62,  the  dipole  moment  p  is  5.62  x  10_^2  c  *cm,  which  is 
equivalent  to  a  charge  separation  distance  of  3.5  A. 

Numerical  solutions  to  the  steady  state  are  presented  below,  afte*-  the 
theoretical  derivation  for  pulsed  operation  of  the  amplifier  is  considered. 

C.  Pulsed  Operation 

For  pulses  with  rise  times  of  the  order  of  1  psec  or  greater,  examination 
of  Eqs.  (16)  through  (22)  and  a  little  algebra  give  a  formal  solution  for  the 
population  difference  N  as 


x  2  x* 

N(k,z\x)  -  {»  +  (Y  N  -  Y8Q0)  f  dx»  exp[(Y  +  Y  )f  +  |^|2dx“]} 

J  0  y2i  Jo 


exP-[(Yp  +  Y8)t  +£j££lrr  |  $  1 2dx 1  ] 


(29) 


Between  pulses,  with  zero  electric  field,  the  solution  for  N  becomes 


n(k,z’,t-t  )  -  (n  (T  )  +y,)(T  -  T  )  -  1]) 

p's  y 


exp  -[ ( Y  +  Y  )(t  -  r  )] 
p  s  o 


where 


N  -  p  -  p 
o  2o  io 


V  P2o  +  plo 


and  Nint  is  the  population  difference  at  the  beginning  of  the  pulse  or  at 
time  At  time  zero  this  population  difference  is 


Nint(k»z'’0)  P2°  "Pl° 


If  the  complex  field  <f'  is  defined  as 


V  m &  exp  i$ 


then  Eqs.  (16)  and  (17)  are  written  as 
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III.  NUMERICAL  RESULTS  AND  DISCUSSION 


A.  Steady  State 

Table  I  shows  the  parameters  used  to  determine  the  semiconductor  device 

U 

properties.  We  have  also  taken  a  typical  spontaneous  emission  rate  y  as  10 

s 

sec-*,  the  equivalent  pumping  rate  as  5  *  10^  sec-*  [17],  and  the  intraband 
relaxation  rate  y^^  as  10*"*  sec-*  [13].  Table  II  gives  the  input  field 
frequencies,  input  field  power  (intensity),  semiconductor  bulk  loss,  and 
assumed  device  series  electrical  resistance  for  each  case  presented  in  this 
report.  Note  that  in  the  results  presented  the  words  intensity  and  power  are 
interchangeable,  because  the  units  of  intensity  are  in  watts  per  square 
micrometer  and  1  pm  of  output  area  is  assumed.  The  line  width  (at  e  )  or 
tne  Gaussian  profile  used  to  cut  off  the  Lorentzian  profile  tails  was  fatten  as 
2  x  10*^  sec-*  for  most  of  the  cases  presented.  Exceptions  are  noted.  As  a 
convenience,  the  calculated  small  signal  gain  for  each  case  is  listed  in 
Table  II. 


The  steady-state  cases  presented  were  calculated  using  a  fourth-order 

Runge-Kutta  algorithm  for  the  differential  Eqs.  (24)  and  (25)  and  by  using  a 

two-point  Gaussian  quadrature  with  automatic  accuracy  checking  to  evaluate  the 

integrals  appearing  in  those  equations.  Rapid  calculation  was  made  possible 

by  evaluating  the  integrals  at  fixed  values  of  dimensionless  intensity  (Eq. 

1/2 

( 26)  J  and  approximating  the  results  by  fourth-degree  polynomial  in  (1+1)  , 

including  negative  powers  of  1/2  and  1. 

Examples  for  steady  state  were  chosen  to  demonstrate  the  effect  of  input 
optical  field  frequency,  semiconductor  bulk  loss,  the  use  of  the  Gaussian 
profile  to  cut  off  the  Lorentzian  line  tails,  and  the  effect  of  inserting  an 
electrical  resistance  in  series  with  the  semiconductor  device.  Figures  1  and 
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LENGTH,  jim 


Steady-state  intensity  response  of  a  GaAs  amplifier.  Input 
parameters  are  listed  in  Table  II.  For  case  1,  use  dBm  scale 
on  right.  For  cases  2,  3,  4  and  5,  scale  factors  are  160,  160, 
15,  and  160  mW,  respectively.  Case  5b  has  1/2-ohm  resistance  in 
series  with  Junction. 


2  show  the  intensity  (power)  and  phase,  respectively,  for  the  five  steady- 
state  cases  listed  in  Table  II.  The  vertical  scales  have  been  normalized  so 
that  all  five  cases  may  be  presented  in  the  same  figure.  In  Fig.  1,  a  decibel 
scale  is  used  for  case  1,  whereas  for  the  other  cases  a  linear  scale  is  used. 
Case  1  is  an  example  of  the  small  signal  regime.  The  frequency  chosen  corre¬ 
sponds  exactly  to  the  band  gap  energy  and  is  at  the  lower  end  of  the  frequency 
range  for  amplification  expressed  by  tq.  (2).  Without  intraband  relaxation 
processes,  the  gain  at  this  frequency  would  be  zero;  however,  collisions  per¬ 
mit  contributions  to  the  gain  at  k  equal  to  zero  from  transitions  at  higher  k 
states.  As  seen  from  Fig.  1,  the  intensity  is  linear  with  device  length,  as 
expected  for  a  decibel  scale,  until  the  intensity  reaches  the  2  to  3 

O 

mW/ pm  range,  at  which  point  saturation  effects  occur.  Saturation  effects  are 
also  apparent  tor  the  phase  shown  in  Fig.  2.  here  the  phase  first  increases 
linearly,  then  the  slope  decreases  as  saturation  occurs. 

We  also  considered  the  other  end  of  the  so-called  allowed  frequency  range 
for  amplification,  at  the  frequency  corresponding  to  the  difference  of  quasi- 
Fermi  levels.  Although  it  is  not  shown  as  a  case,  the  gain  was  calculated  to 
be  negative.  If  tq.  (2)  were  strictly  valid,  this  gain  would  be  zero;  how¬ 
ever,  because  of  a  finite  pumping  rate  and  intraband  relaxation,  a  fairly 
large  absorption  results.  Case  2  is  calculated  at  a  transition  frequency 
giving  the  largest  small  signal  gain.  For  this  case,  Fig.  1  shows  a  rapid 
transition  from  the  small  signal  regime  to  a  linear  gain  region.  The  phase, 
shown  in  Fig.  2,  is  negative  for  this  case,  which  is  not  surprising  because  we 
have  convolved  the  population  difference  with  the  standard  expression  for  the 
real  part  of  the  susceptibility  (Ref.  28,  p.  154)  of  a  single  state.  This 


LENGTH,  pm 


Steady-state  phase  response  of  a  GaAs  amplifier.  Scale  factors 
for  cases  lf  3,  4,  and  5  are  14,  0.17,  5,  and  5  rad,  respectively 
for  case  2,  scale  factor  on  the  right  Is  4  rad. 


susceptibility  is  known  to  be  antisymmetric  about  the  resonance.  If  there 
were  no  upper  limit  on  the  k  integration,  or  if  the  Gaussian  cutoff  were  not 
used,  this  phase  would  approach  negative  Infinity.  A  switch  in  sign  tor  the 
phase  between  cases  1  and  2  indicates  that  there  is  a  frequency  that  yields  a 
zero  phase  change  for  the  amplifier.  lhat  this  is  true  is  shown  for  case  3, 
with  an  input  frequency  of  A.  28  x  10^  hz.  The  intensity  for  case  3  shown  in 
Fig.  1  is  similar  to  that  for  case  2  but  with  a  slightly  smaller  gain.  The 
phase  is  very  interesting  in  that  there  is  a  sign  change  with  device  length. 
As  the  intensity  of  the  optical  field  grows  in  the  amplifier,  the  effective 
homogeneous  line  width  also  increases.  This  permits  contributions  of  higher  k 
states  resulting  in  a  net  decrease  in  the  phase  rather  than  an  increase  as  in 
the  small  signal  regime.  Note  that  the  frequencies  for  maximum  gain  and  zero 
phase  change  are  not  the  same. 

Case  4  demonstrates  strong  saturation  effects.  The  frequency  is  close 
but  not  exactly  equal  to  that  used  for  the  examples  for  the  pulsed  operation 
of  the  amplifier.  This  frequency  was  chosen  to  be  slightly  different  than  for 
the  zero  phase  case,  because  zero  phase  requires  longer  computation  times. 
The  hulk  loss  was  increased  to  show  saturation.  The  intensity  shown  in  Fig.  1 
saturates  rapidly  to  a  value  of  14  mW/jm^.  The  phase  for  this  case  exhibits 
two  linear  regimes;  the  first  is  for  the  small  gain  regime  where  the  intensity 
is  close  to  zero,  and  the  second  is  in  the  region  where  the  intensity  has 
saturated  to  its  final  value. 

For  case  5,  we  considered  the  effect  of  placing  an  electrical  resistance 
in  series  with  the  semiconductor  device  to  limit  the  current  drawn  by  the 


semiconductor  during  amplification.  This  resistance  may  either  he  intrinsic 
to  the  device  or  included  in  the  associated  external  circuit.  The  bias 
voltage  across  the  junction  drops  during  amplification  because  of  this  series 
resistance,  permitting  a  decrease  in  the  quasi-Fermi  levels,  which  implies 
that  the  device  will  exhibit  less  gain  than  for  the  unrestricted  or  no 
resistance  case.  The  effect  of  placing  a  1/2-ohm  resistance  in  series  with 
the  aiq>llfler  is  shown  as  case  5b  in  Figs.  1  and  2.  The  corresponding  no 
resistance  case  is  5a.  As  expected,  the  small  signal  gain  decreases,  but  only 
slightly.  On  the  other  hand,  the  phase  decreases  by  a  factor  of  2,  and  it 
appears  that  the  phase  saturates,  whereas  for  the  noncurrent-limited  case 
phase  continues  to  increase  with  the  length  of  the  amplifier. 

To  test  the  model  dependence  on  the  various  parameters  we  have  used, 

several  cases  were  run  with  parameters  different  from  those  given  above. 

These  cases  are  shown  in  Figs.  3  and  A  and  the  parameters  used  are  indicated 

in  Table  III.  Cases  presented  are  labeled  1,  2a,  2b,  and  2c.  Cases  1  and  2a 

are  to  be  compared  with  cases  1  and  2  above,  respectively,  whereas  cases  2b 

and  2c  are  to  be  compared  with  each  other,  case  2a,  and  case  2  above.  Note 

that  cases  2a,  b,  and  c  are  for  the  same  frequency.  Cases  1  and  2a  were 

calculated  with  a  collisional  relaxation  rate  2.5  times  higher  than  used 

previously.  Corresponding  to  this  relaxation  rate,  the  Gaussian  linewidth  was 

14 

taken  as  20  Tr21  or  5  x  10  Hz.  The  increase  in  linewidth  for  both  the 
Lorentzlan  and  Gaussian  profiles  permits  a  larger  range  of  k-states  to  contri¬ 
bute  to  the  gain  and  phase  Integrals  lEqs.  (24)  and  (25)].  This  means  that 
the  small  signal  gain  will  increase  for  case  1  and  decrease  for  case  2a 
compared  to  that  seen  for  cases  1  and  2  above,  respectively.  This  is  verified 
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POWER,  mW 


LENGTH,  im 

Steady-state  Intensity  response  of  a  GaAs  amplifier.  For  case  1, 
use  dBm  scale  on  right.  For  cases  2,  3,  and  4,  scale  factors  are 
70,  70,  and  240  nW,  respectively.  Table  III  shows  input  values. 


Table  III 

ALTERNATE  OPTICAL  INPUT  PARAMETERS  AND  CALCULATED  SMALL  SIGNAL  GAIN 


Steady  state 


by  calculating  the  small  signal  gains  listed  in  Table  111.  The  gain  for  case 
2a  decreases,  because  the  inclusion  of  more  k-states  in  the  gain  integral 
averages  in  states  of  lower  gain.  The  inclusion  of  lower  gain  k-states 
Increases  with  Intensity,  resulting  in  intensity  saturation  at  a  value 
of  70  mW/um  ,  as  can  be  seen  from  Fig.  3.  Phase  is  also  affected  by  the 
Inclusion  of  added  k-states.  The  phase  for  case  1  reaches  a  maximum  and 
decreases,  which  is  similar  to  the  phase  behavior  for  case  3,  although  the 
peak  phase  is  considerably  larger  (3.5  rather  than  0.17  radians).  The  phase 
for  case  2a  decreases  much  more  rapidly  than  for  case  2,  reaching  a  value  of 
nearly  -30  radians  as  compared  to  -4  radians  (Fig.  4). 

To  verify  a  prediction  of  very  large  phase  in  the  absence  of  the  Gaussian 
cutoff,  we  recalculated  case  2  (case  2b)  above  for  a  relaxation  rate  of  10*"* 
Hz  but  eliminated  the  Gaussian  cutoff.  As  expected  the  phase  decreases 

rapidly— 12  times  faster  than  previously.  The  intensity  also  shows  saturation 

effects  caused  by  the  removal  of  this  cutoff.  The  last  case  shown,  2c,  was  a 

modification  of  2b.  The  relaxation  rate  was  increased  by  a  factor  of  10  to 

10U  Hz  but  the  range  of  k-states  permitted  to  contribute  to  the  gain  and 
phase  was  restricted  by  lowering  the  upper  limit  on  the  k  integration  by  30% 
to  6.47  x  10^  cm-*.  This  k  cutoff  corresponds  to  that  value  of  k  for  which 
the  direct  (D  and  indirect  (X)  band  energies  are  the  same.  The  small  signal 
gain  for  2c  decreases  with  respect  to  2b  by  about  30%.  Note  that  the  scale 
factor  for  phase  in  Fig.  4  is  only  0.3  for  case  2c,  indicating  a  very  small 
phase  variation  for  this  case  because  of  the  nearly  symmetric  inclusion  oc 
k-states  with  positive  and  negative  contributions  to  the  phase  Integral 
(bq.  25).  Case  5  (not  shown)  with  and  without  a  series  resistance  was 
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recalculated  for  *  10  /sec,  no  Gaussian  tail  cutoff,  and  the  k  upper 

limit  of  case  4.  The  significant  behavior  difference  was  a  decreased  small 
signal  gain  and  a  phase  change  from  positive  to  negative  without  and  with  the 
series  resistance,  respectively.  A  shift  in  population  level  with  bias 
voltage  explains  the  variation. 

The  model  dependence  on  the  chosen  parameters  indicates  that  experiments 
snould  be  made  to  verify  results  presented  here, 
d.  Pulsed  Operation 

Examples  presented  here  demonstrate  the  effect  of  the  amplifier  on  a 
pulse  or  a  series  of  pulses.  Pulse  compression,  phase  nonuniformity  with 
respect  to  time  for  the  output,  and  multi-pulse  effects  are  of  interest.  The 
assumed  pulse  shape  is  shown  as  an  insert  in  Fig.  5.  The  input  pulses  were 
chosen  to  be  54  psec  long.  A  sinusoid  electric  field  with  a  3  x  1010  v/m  peak 
(peak  intensity  0.44  mW/pm  )  at  29  psec  and  a  half  cycle  length  of  50  psec 
were  used  as  the  basic  input.  The  basic  shape  was  modified  by  providing  a 
linear  ramp  from  zero  time  such  that  this  ramp  was  tangent  to  the  sine  at 
approximately  17  psec.  The  Gaussian  shape  is  not  causal  and  was  not  con¬ 
sidered.  Note  that  the  exact  pulse  shape  is  not  important  for  determining 
qualitative  features  of  the  model.  Optical  input  parameters  used  for  the 
calculation  are  shown  in  Table  II  as  cases  6  and  7. 

The  numerical  procedures  used  to  solve  Eq.  (35)  differ  considerably  from 
the  techniques  used  for  the  steady-state  case.  A  high-order  Runge-Kutta 
algorithm  was  used  to  perform  the  z  integration,  whereas  the  trapezoidal  rule 
was  used  to  perform  the  time  and  momentum  integration.  The  steps  in  the 
calculation  follow.  At  each  z  step  the  electromagnetic  field  as  a  function  of 
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Fig.  5.  GaAs  amplifier  Intensity  response  to  a  single  input  modified 

sinusoidal  pulse  for  various  device  lengths.  The  Insert  shows  the 
input  pulse  amplitude  shape. 


time  is  given  by  the  results  at  the  previous  value  of  z.  From  these  values, 
the  population  inversion  as  a  function  of  time  and  k  is  determined  from  Eqs. 
(29)  and  (30).  The  required  input  to  the  Runge-Kutta  algorithm  is  completed 
by  performing  the  momentum  integration  in  Eq.  (35). 

Pulse  compression  and  output  phase  variations  with  time  are  shown  in 
Figs.  5  and  6.  The  frequency  was  chosen  to  obtain  a  high  gain  but  a  rela¬ 
tively  small  phase  change  during  amplification.  In  these  figures  the  curves 
presented  are  Intensity  or  phase  versus  time  for  various  values  of  device 
length  as  a  parameter.  The  pulse  compression  shown  in  Fig.  5  is  extreme 
because  of  the  high  small  signal  gain.  The  leading  edge  of  the  pulse  expe¬ 
riences  a  gain  of  a  factor  of  55  for  each  50-um  length,  which  results  in  a 
pulse  compression  to  2  psec  in  an  amplifier  length  of  200  um.  Compression 
occurs  because  the  population  density  is  a  function  of  the  integrated  inten¬ 
sity  [Eq.  (29)],  so  that  a  high  intensity  results  in  reduced  gain  or  net 
absorption  for  the  remaining  pulse.  Compressions  shown  that  are  less  than 
2  psec  are  probably  not  valid  because  approximation  9  no  longer  holds  and  the 
complete  set  of  Eqs.  (3)  through  (6)  and  (12)  must  be  solved  [10],  Note 
that  we  have  attempted  to  limit  the  phase  variation  with  time  by  a  choice  of 
frequency.  Relative  success  can  be  seen  in  Fig.  6,  which  exhibits  a  phase 
change  of  approximately  0.13  radian  for  the  200-ym  amplifier,  which  still 
amounts  to  a  change  of  30%.  Other  choices  for  frequency  yield  a  phase  change 
of  several  radians. 

In  Section  III-A  we  considered  the  effect  of  a  resistance  in  series  with 
the  heterostructure  junction  for  a  steady-state  case.  The  effect  of  a  series 
resistance  on  pulsed  outputs  is  difficult  to  assess  quantitatively  with  the 
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model  provided  here.  Qualitatively,  amplification  results  in  lowering  the 
bias  across  the  junction.  Gain  is  decreased  while  the  pulse  travels  through 
the  amplifier  resulting  in  less  pulse  compression  and  a  more  uniform  phase 
than  shown  for  Figs.  5  through  8.  It  would  be  interesting  to  model  and 
calculate  the  effect  of  such  a  resistance  as  well  as  to  perform  experiments 
verifying  such  predictions. 

C.  Multi-Pulse  Effects 

For  applications  that  require  amplification  of  a  sequence  of  pulses  such 

as  occur  in  an  optical  communication  system,  knowledge  of  multiple-pulse 

effects  is  of  primary  interest.  For  this  study,  six  identical  pulses  were 

amplified  in  sequence  by  a  200-um  long  device;  the  resulting  multipulse 

effects  are  exhibited  in  Figs.  7  and  8.  A  precursor  steady-state  intensity 
2 

of  0.01  mW/um  was  input  to  reduce  the  population  inversion  before  pulse 
arrival.  This  precursor  was  turned  off  50  psec  before  the  first  pulse 
arrived.  The  pulse  sequence  is  illustrated  in  Fig.  7.  The  differences  in 
output  for  the  various  pulses  in  the  sequence  are  caused  by  the  slow  pumping 
time  in  relation  to  the  pulse  length  and  interpulse  spacing.  Note  that  there 
is  significant  inversion  recovery  between  the  third  and  fourth  pulses. 
Considerable  pulse-to-pulse  variation  occurs  in  phase  as  well  as  in  inten¬ 
sity.  For  signal  processing  applications,  this  phase  variation  oust  be 
minimized  by  frequency  choice. 

From  assumption  9,  pulse-to-pulse  interference,  per  se,  is  not  possible 
when  the  pulses  do  not  overlap.  Therefore,  after  a  pulse  has  passed  through 
the  amplifier,  there  is  no  residual  polarization  that  can  generate  an  optical 
field  (ringing  or  photon  echoes)  as  is  possible  when  the  full  set  of  equations 
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GaAs  amplifier  phase  response  to  a  set  of  six  input  pulses 
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is  considered.  Therefore  each  input  pulse  generates  an  individual  output 
pulse  based  entirely  on  the  inverted  population  levels  at  the  beginning  of  the 
pulse  and  the  input  pulse  shape.  However,  because  the  pumping  processes 
generating  the  inverted  population  levels  are  much  slower  than  the  collisional 
processes,  the  GaAs  amplifier  exhibits  apparent  interpulse  interference  for 
repetition  rates  approaching  10  GHz,  as  shown  in  the  figures.  This  occurs 
because  there  is  insufficient  time  between  pulses  to  build  up  to  the  signal- 
free  Inverted  population  levels.  These  effects  can  be  avoided  by  limiting 
saturation  in  the  amplifier  by  shortening  the  amplifier  length  or  decreasing 
the  available  gain. 

From  the  analysis  presented,  devices  similar  to  those  modeled  here  would 
amplify,  compress,  and  show  a  phase  variation  for  every  pulse  at  a  1  Gbit  rate 
for  50  psec  pulses.  These  pulses  would  be  similar  to  those  shown  in  Figs. 
5  and  6.  Amplification  of  a  continuous  pulse  train  at  10  Gbits  would  show 
results  similar  to  those  shown  for  cases  1,  2,  5,  and  6  of  Figs.  7  and  8 
unless  at  least  three  pulses  were  skipped. 

For  low  data  rates,  phenomena  such  as  pulse  compression  and  phase  change 
during  the  pulse  must  be  compensated  in  applications  that  require  uniform 
intensity  or  phase.  For  high  data  rates,  a  constant  bit  stream  will  keep  the 
Inverted  population  at  a  low  level  and  produce  a  uniform  output.  For  communi¬ 
cations  applications,  a  constant  bit  stream  could  encompass  coding  schemes 
such  as  polarization  modulation  or  pulse  position  modulation  in  which  the 
location  of  a  pulse  is  shifted  in  time  only  a  short  distance  but  never 
deleted.  On-off  coding,  on  the  other  hand,  skips  pulses,  thereby  permitting 


the  Inverted  population  to  grow,  leading  to  nonuniform  shapes  pulse-to-pulse 

therefore,  we  suggest  the  use  of  the  former  coding  schemes  when  pulse  uni 
fortuity  is  required. 


IV.  CONCLUSIONS 


Amplification  in  monolithic  GaAs  by  means  of  a  self-consistent  semi- 
classical  approach  has  been  considered.  Phenomenological  rate  constants  were 
Included  to  account  for  spontaneous  emission,  pumping,  and  collision-induced 
decay.  Their  magnitude  was  based  on  experimental  evidence  concerning  the 
homogeneous  line  width  and  electron-electron  induced  decay  times.  The  pumping 
rate  constant  was  chosen  to  provide  strong  evidence  of  saturation  effects  in 
steady  state  in  the  region  of  10  mU  output,  which  seems  to  be  justified 
experimentally  130]. 

The  results  obtained  must  be  considered  qualitative  and  not  quantitative, 
because  experiment  has  provided  neither  definitive  values  for  the  collisional 
relaxation  rates  nor  the  exact  shape  of  the  collision-induced  emission  line 
profile.  These  parameters  must  be  determined  before  an  accurate  calculation 
of  gain  and  phase  change  through  an  amplifier  can  be  performed.  Since  the 
results  presented  here  are  model  dependent,  experimental  results  are  needed  to 
predict  the  behavior  of  the  phase  and  gain  in  amplifiers. 

The  qualitative  information  available  from  the  analysis  is  extensive. 
For  the  assumed  range  of  collisional  relaxation  rates  (10l3/sec  to  lO^/sec), 
the  use  of  detailed  rate  equations  to  describe  optical  amplification  was  found 
to  be  appropriate  for  data  rates  as  high  as  10  GHz  and  input  pulse  widths  as 
short  as  20  psec,  as  long  as  the  pulses  do  not  compress  to  less  than  2  psec 
during  amplification.  Along  with  a  rate  equation  for  Intensity,  we  have  given 
a  new  rate  equation  for  phase,  which  can  be  used  in  parallel  with  the  equation 
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for  gain.  These  rate  equations  are  valid,  with  the  assumption  that  the  polar 
ization  follows  the  electric  field  exactly. 
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APPENDIX 


NONREGENERATIVE  OPTICAL  AMPLIFICATION 
IN  MONOLITHIC  GaAs 


1.  INTRODUCTION 

In  this  Appendix,  a  theoretical  investigation  of  nonregeneratlve  light 
amplification  in  GaAs  double  heterostructure  devices  is  presented.  The  analy¬ 
sis  is  carried  out  from  first  principles.  Previous  analyses  of  GaAs  devices 
were  based  on  simple  rate  equations  for  laser  devices.  These  rate  equations 
do  not  contain  information  on  the  phase  behavior  of  pulses  propagated  through 
the  device  nor  on  the  range  of  data  rates  for  which  the  theory  is  valid. 

This  investigation  was  carried  out  in  conjunction  with  an  experimental 
program  that  had  the  objective  of  developing  a  GaAs  amplifier  capable  of  a 
coherent,  small  signal  power  gain  of  1000  and  coherent  large  signal  behavior. 
Further  goals  were  high  data  rate  capability  without  pulse-to-pulse  interfer¬ 
ence  effects  and  a  large  frequency  range  for  amplification.  The  final  goal  of 
a  large  frequency  bandwidth  is  not  practical  for  laser  amplifiers  because  of 
the  Fabry-Perot  or  distributed  Bragg  reflector  conditions  that  determine  the 
frequency  of  operation.  For  these  reasons,  it  was  decided  to  examine  nonre- 
generative  optical  amplif ication  in  GaAs  from  first  principles.  Nonregenera- 
tive  devices  are  fabricated  to  prevent  optical  feedback  and  as  such  do  not 
self-osclllate.  Since  there  are  no  oscillation  conditions,  any  frequency 
within  the  spontaneous  emission  bandwidth  of  GaAs  can  be  amplified.  Further¬ 
more,  plane  wave,  rather  than  modal  analysis,  is  used  since  the  optical  wave 
passes  through  the  device  only  once. 

In  Section  2,  a  detailed  exposition  of  the  theory  begins,  after  a  brief 
introduction,  with  the  reduction  of  the  material  characteristics  to  an  equiva¬ 
lent  two-level  atom  problem,  the  introduction  of  the  density  matrix  formation, 
and  a  discussion  of  the  phenomenologically  introduced  decay  processes.  A 
brief  discussion  of  the  homogeneous  and  inhomogeneous  line  widths  is  also 
given.  In  Section  3,  simplifying  assumptions  are  made,  the  interaction 


Hamiltonian  is  introduced  and  quantized,  and  the  equations  of  motion  for  the 
density  material  and  electromagnetic  field  are  derived.  Finally,  the  equa¬ 
tions  are  simplified  and  reduced  to  rate  equations  because  of  the  large  value 
for  the  collision-induced  decay  processes.  In  Sections  4  and  5,  the  equations 
are  specialized  to  steady-state  or  pulsed  conditions,  respectively.  Device 
parameters  and  an  explanation  of  the  derivation  leading  to  the  chosen  device 
parameters  are  given  in  Section  6.  The  specialization  to  a  given  set  of 
device  parameters  is  necessary  to  carry  out  numerical  results,  which  are 
Included  in  the  main  body  of  the  text. 

2.  THEORY 

The  investigation  presented  in  this  appendix  is  based  on  a  self- 
consistent,  8emidassical  approach.  The  electromagnetic  field  is  treated 
classically,  while  the  medium  is  quantized.  Self-consistency  is  Invoked  by 
requiring  that  the  electric  field  Introduce  a  polarization  in  the  medium 
which,  in  turn,  acts  as  a  source  term  to  Maxwell's  equation  (Fig.  A-l).  The 
electromagnetic  field  is  assumed  to  be  a  linearly  polarized  plane  wave.  The 
medium  consists  of  a  group  of  two-level  electron  states  generated  by  applying 
a  bias  voltage  across  a  monolithic  P-n-N  he teroj unction  GaAs  device.  The 
conduction  and  valence  band  occupations  are  described  by  Fermi-Dirac  distri¬ 
bution  functions.  Electron  transitions  between  the  two  energy-level  dis¬ 
tributions  rather  than  between  two  discrete  levels  must  be  considered  for 
interband  absorption  or  emission  of  radiation  at  a  given  frequency.  The 
dipole  approximation  is  used  to  describe  the  interaction  between  Che  electric 
field  and  two  level  states.  Phenomenological  rate  constants  are  Included  in 
the  equations  to  describe  the  effect  of  pumping,  spontaneous  emission,  and 
photon  or  electron-electron  collision-induced  decay.  The  effect  of  sponta¬ 
neous  emission  noise  is  not  Included  in  the  analysis  although  an  understanding 
of  this  phenomenon  is  necessary  for  a  complete  solution  to  the  problem.  In 
order  to  quantify  this  problem  somewhat,  an  actual  device  is  contemplated. 


Let  the  device  consist  of  three  layers,  P-n-N  (shown  in  Fig.  A-2),  with  doping 
and  aluminum  concentrations  as  indicated.  The  unbiased  device  is  shown  in 

Fig.  A-2(A),  and  the  biased  device  is  shown  in  Fig.  A-2(B).  These  figures  are 
not  drawn  to  scale.  The  details  leading  to  these  parameters  are  shown  in 
Section  6  and  are  taken  from  information  obtained  in  Ref.  A-l.  For  simplicity 
it  is  assumed  that  there  are  no  "band  tails"  and  that  the  k-selection  rules 
hold.  This  establishes  the  relationship  between  the  circular  frequency  w  of  a 
given  transition  and  its  corresponding  wavenumber  k  as 

■fc2k2 

*  Eg  +  2bT~  (A-1) 

where  Eg  is  the  band-gap  energy  (1.6734  eV)  of  intrinsic  Ga80Al20*®  at  room 
temperature  shown  in  Fig.  A-2,  mr  jg  the  reduced  ness  (7.2428  x  10“  2  m0  = 

6.596  x  10"29  g)t  and  is  Planck's  constant  divided  by  2ir.  A  schematic 
representation  of  the  direct  bands,  assuming  a  parabolic  dependence  on  k,  is 
shown  in  Fig.  A-3.  Note  that  this  schematic  does  not  correspond  exactly  to 
that  shown  in  Fig.  A-2,  but  corresponds  more  to  the  classic  textbook  repre¬ 
sentation  of  band  models.  The  purpose  of  Fig.  A-3  is  to  point  out  the 
parabolic  nature  of  the  bands  and  indicate  the  shape  of  the  Fermi  distribu¬ 
tions. 

The  theoretical  treatment  follows  that  presented  by  A.  Icsevgl  and  W.  E. 
Lamb,  Jr.  (Ref.  A-2)  for  a  discrete  atomic  system.  Since  there  are  differ¬ 
ences  between  problems  with  discrete  atomic  transitions  and  transitions 
between  energy  bands,  we  develop  the  theory  from  first  principles. 

As  stated  above,  the  medium  is  considered  to  be  modeled  using  two  energy 
bands  (Fig.  A-3)  coupled  by  their  interaction  with  an  overall  radiation 
field.  If  a  population  inversion  between  the  conduction  band  and  valence  band 
can  be  established,  the  medium  is  capable  of  amplifying  light  in  a  frequency 
band  determined  by  the  band-gap  and  quasi-Fermi  levels. 

E  <  liw  <  E„  -  E  (A-2) 
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A-2.  P-n-N  Heterojunction  Device  Without  and  With  Bias  Voltage 


In  order  to  carry  out  necessary  statistical  summations  to  eliminate 
unobserved  variables,  it  is  convenient  to  represent  the  conduction  and  valence 
bands  as  if  they  were  levels  of  a  two-level  atom,  thus  obtaining  equations  for 
the  density  of  states  in  either  band.  Following  the  notation  for  a  two-level 
atom,  it  is  convenient  to  introduce  a  2  by  2  density  matrix  p.  The  diagonal 
elements  will  be  related  to  the  density  of  states  in  the  conduction  and 
valence  bands.  The  upper  state  represents  an  electron  in  the  conduction  band 
and  a  corresponding  hole  in  the  valence  band,  whereas  the  lower  state  repre¬ 
sents  an  electron  in  the  valence  band  and  a  corresponding  hole  in  the  con¬ 
duction  band.  The  off-diagonal  elements  are  related  to  the  dipole  moment. 

The  density  matrix  is 


The  decay  of  levels  c  and  v,  as  well  as  the  damping  of  the  dipole  moment,  must 
be  described  phenomenologically  by  introducing  such  constants  as  Yg,  Yco^, 
and  Ycv«  where  ycQ^  represents  the  decay  of  levels  resulting  from  collision- 
induced  decay,  the  quantity  y  corresponds  to  the  spontaneous  emission  rate, 

S 

and  Ycv  corresponds  to  dephasing. 

Induced  decay  causes  a  population  level  ”n"  to  decay  as  follows: 

n(t)  ■  nj^  exp(-  y  jt)  (A-4) 

where  n^  is  the  initial  population  density.  At  the  same  time,  pumping  will 
establish  the  population  inversion.  The  spontaneous  emission  time  is  of  the 
order  of  10-9  sec,  while  the  collision  damping  constant  is  of  the  order  of  the 
inverse  electron-electron  collision  time.  Since  the  "phase**  of  the  dipole 
moment  of  a  transition  is  lost  when  one  of  the  upper  or  lower  states  is 
destroyed  as  a  result  of  collision,  one  could  say  that  the  decay  time  is  es¬ 
sentially  given  by  Associated  with  ycv  is  the  (homogeneous)  line  width 

of  the  emission.  An  additional  broadening  exists  because  of  the  band  nature 
of  the  conduction  and  valence  bands.  This  establishes  the  inhomogeneous  line 


width.  The  entire  width  of  the  fluorescence  emission  is  given  by  a  convo¬ 
lution  of  the  homogeneous  and  inhomogeneous  line  width.  Investigations  of 
electron-electron  collision  times  (Ref.  A-3)  and  the  homogeneous  line  width 
(Ref.  A-4),  as  well  as  private  discussion,*  have  Indicated  a  value  of  and 

therefore  y  of  10^  to  10^.  For  calculations  presented  herein  we  have 
CV  13 

taken  y  ,  *  Y  *  10  .  In  the  next  section,  the  details  of  the  theoretical 

col  cv  * 

investigation  are  considered.  Phenomenological  terms  corresponding  to  pumping 
are  also  Included  there. 

3.  THEORETICAL  DEVELOPMENT 

A.  ASSUMPTIONS 

The  calculation  to  be  performed  is  summarized  in  Fig.  A-l.  The  assumed 
electric  field  £(r,t)  polarizes  the  electron  in  the  medium.  The  electronic 
dipole  moments  are  summed  to  give  the  macroscopic  polarization  $(r,t),  which 
enters  Maxwell's  equation  as  a  source  term.  The  electric  field  driven  by  this 
source  term  must  then  equal  the  applied  electric  field  to  be  self-consistent. 
One  begins  by  determining  the  interaction  Hamiltonian  for  the  two-level 
states.  The  two-level  system  is  quantized,  and  the  equations  of  motion  for 
the  density  matrix  are  determined.  Taking  into  account  the  band  nature  of  the 
amplifying  medium,  the  expressions  are  summed  to  give  the  macroscopic  polari¬ 
zation.  Finally  this  polarization  is  introduced  back  into  Maxwell's  equa¬ 
tions.  In  order  to  derive  these  equations,  the  following  simplifications  will 
be  applied: 

1.  The  two-level  states  are  coupled  only  through  their  interaction  with 
the  common  field  £(r,t). 

2.  The  radiation  field  is  a  scalar  E(z,t)  representing  a  uniform  plane 
wave  polarized  in  only  one  direction.  We  take  this  polarization 
direction  as  x  with  propagation  in  the  z  direction.  This  means  the 
questions  concerning  regeneration  and  modal  structure  are  completely 
Ignored. 

3.  The  dipole  approximation  holds  for  the  interaction  of  the  electrons 
with  the  radiation  field. 


Haus,  private  communication. 
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4. 


Amplification  in  a  double  heterostructure  device  is  being 
considered.  The  lasing  region  is  assumed  to  be  intrinsic  GaAlAs. 

For  this  reason  the  parabolic  band  approximation  is  used  and  the 
only  allowed  transitions  obey  the  k  selection  rule,  i.e. , 
transitions  are  between  states  with  exactly  the  same  momentum. 

B.  INTERACTION  HAMILTONIAN 

The  classical  form  of  the  nonrelativistic  Hamiltonian  for  a  single 
particle  of  mass  m  and  charge  q  moving  in  an  electromagnetic  field  is  given  by 

H  -  jg  (p  -  q£)2  +  q  ♦  (A- 5) 


where  p  is  the  momentum  conjugate  to  the  position  of  the  particle,  A  is  the 

'0' 

magnetic  vector  potential,  and  4>  is  the  scalar  field.  The  Hamiltonian  is 
derived  in  the  usual  way  (Ref.  A-5),  i.e.,  from  a  classical  Lagrangian  that 
includes  a  generalized  velocity-dependent  potential  obtained  from  the  Lorentz 
force  on  the  particle  arising  from  the  electromagnetic  field  (1,$).  That 
field  obeys  Maxwell's  equations,  given  here  for  completeness: 


(A-6) 

5  -  p  -  o 

(A-7 ) 

(A-8) 

•  I  -  0 

(A-9) 

Equation  (A-9)  implies  that 


and  Eq.  (A-6)  that 


S-^xX 


(A-10) 
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C A— 1 1) 


The  conversion  from  classical  physics  Co  quantum  physics  is  accomplished  by 
the  replacement  of  the  conjugate  momentum  p  by  ti/i$.  For  convenience  p  will 
be  retained  wherever  possible.  An  additional  potential  energy  V  is  added  to 
the  Hamiltonian  in  Eq.  (A-5). 

The  single-particle-*plus-fleld  Hamiltonian  is  written  in  the  form 


H  -  Hm  +  Hc  +  Hi 


CA-12; 


where 


“„--k?2  +  v 


(A-13) 


HI  "  m  1  '  P  +l^  *  X)  +  2S  CAl2  +  q* 


(A-14) 


and  the  final  term  is  the  energy  content  of  the  electromagnetic  field 


Hf  -  j  /  (S  •  $  +  $  •  !)d3x 


(A-15) 


Since  the  electromagnetic  strength  does  not  vary  over  distances  of  the  order  . 
of  a  wavelength,  it  is  possible  to  sum  the  interaction  energy  over  the  allowed 
states  in  the  conduction  band  with  the  states  labeled  by  the  momentum  k 


Hx(x) 


I  (i  v  -h+T^v-  V) 


(A-16) 


By  performing  a  gauge  transformation,  one  obtains 


(A- 17) 


The  quantity  qx^  is  recognized  as  the  dipole  moment  therefore 

Hj  -  2(z)  •  l  wk 


(A-24) 


C.  EQUATIONS  OF  MOTION  FOR  THE  DENSITY  MATRIX 

It  is  well  known  that  the  equation  of  motion  for  the  density  matrix  for  a 
single  "pure"  state  labeled  by  k  is  given  by  (Ref.  A-6) 
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* 


(A-32) 


p  “  P 
cv  VC 

where  the  phenomenological  terms  y  ,  y  , ,  and  y  have  been  introduced  to 

8  co  1  cv 

account  for  both  damping  and  pumping.  Note  that  pumping  is  established  by 
applying  a  bias  across  the  GaAs  device.  This  bias  establishes  equilibrium 
values  of  the  upper  and  lower  states  p  and  p  that  would  exist  in  the 
absence  of  spontaneous  emission  and  an  electromagnetic  field.  Differences 
from  equilibrium  then  accomplish  pumping  or  decay  through  the  collision 
process.  These  equations  will  be  modified  slightly  later  in  the  report  but 
are  now  sufficient  to  continue  our  discussion. 

The  individual  electron  or  Bloch  states  may  be  specified  by  various 
labels,  such  as  position  z,  time  t,  momentum  k,  upper  or  lower  state  (c,v), 
etc.  Note  that,  since  we  use  the  envelope  or  slowly  varying  wave  approxi¬ 
mation,  we  may  specify  z  and  k  at  the  same  time.  The  dipole  moment  p  is  given 
by  the  trace  with  the  density,  resulting  in 

<p>  -  p[pcv(z,t,r,  etc.)  +  Pvc]  (A-33) 

On  the  other  hand,  Eqs.  (A-29)~(A-31)  depend  only  on  z,  t,  and  k,  and  the 
macroscopic  polarization  of  the  medium  P(z,t)  does  not  depend  on  k.  For  this 
reason  the  extra  variables  labeling  the  individual  electron  states  are 
eliminated  by  taking  the  ensemble  average  for  all  states  having  the  same  t,  z, 
and  k  labels,  and  the  total  polarization  P(z,t)  is  then  the  sum  over  partial 
polarizations. 

P(z,t)  -  /  dkP(k,z,t)  ■  /  dkp[pcv(k,z,t)  +  cc]  (A-34) 

Equation  (A- 29)  is  rewritten  with  specific  labels  as 
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while  Eqs.  (A-30)  and  (A-31)  are  written  for  completeness  as 


9p(k, z ,  t)  -  ~Y  .  (p  -  P  )  +  Y  P  +  i(-^r)  (P^  -  P„J  (A-36) 

V  col  V  VO  8  C  v  J  cv  VC 

3pcv(k*z‘t)  "  "(YCV  +  iu)pcv  “  1(jr)E<*»t><pc  ■  PV}  (A'37) 

where  9  represents  9/9t. 

These  equations  can  be  written  in  terms  of  the  partial  polarization 
P(k,z,t)  in  a  straightforward  manner  (Ref.  A-2): 

P2 

[O  +  Ycv)2  +  w2]p(k,z,t)  -  -2w(2^)E(z,t)(pc  -  Py)  (A-38) 


O  +  Y  -.  )p  -  Y  |P  +  Y.P„  -  ‘ho)"1  E(9  +  Y  )P  (A-40) 

col  v  col  VO  s  c  cv 

The  corresponding  electromagnetic  equation  follows  immediately  from 
Maxwell' 8  equations,  giving 


32E  ,  a3E  .  2  -2  92E 

~  ~T  +  Vo  “Tt  +  noC  ~2 


Pq3  P(z,t) 


(A-41) 


where  the  conductivity  o  accounts  for  any  linear  losses  incurred  in  the 
absorbing  medium,  c  is  the  velocity  of  light  in  vacuum,  and  is  the  magnetic 
permeability.  Note  that  the  total  polarization  enters  Eq.  (A-41),  while  the 
partial  polarizations  enter  Eqs.  (A-38)-(A-40) »  The  intrinsic  polarization  of 
the  material  resulting  from  the  background  crystal  matrix  is  included  as  a 
factor  n0,  the  index  of  refraction. 
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SIMPLIFICATION  OF  THE  BASIC  EQUATIONS 


We  now  make  two  additional  assumptions  and  remind  the  reader  of  a  third 
already  discussed. 


1.  The  optical  frequency  uj  is  larger  than  the  optical  linewidth; 
therefore 


ui  >  y 


cv 


(A-42) 


2.  The  rotating  wave  approximation  applies,  i.e. ,  terms  involving  the 
creating  of  a  photon  and  the  simultaneous  creation  of  the  upper 
state  leading  to  harmonics  of  the  optical  field  are  ignored. 

3.  Recall  that  the  slowly  varying  wave  approximation  has  already  been 
invoked  in  the  description  of  the  medium.  For  the  electromagnetic 
wave  this  titfans  we  can  write 


E(z,t)  «<£(z,t)  cos[vt  -  nz  +  $(z,t)] 


(A-43) 


The  frequency  circular  v  will  have  a  nominal  value  close  to  the  center  of  the 
bandwidth  of  the  input  pulse.  Different  choices  of  v  lead  to  different  phase 
functions  4>(z,t).  The  wavenumber  n  will  always  be 


n  v 
o 


(A-44) 


The  amplitude  and  $  phase  will  vary  slowly  in  an  optical  cycle  and  wavelength. 


i.e.  , 


If  «  If  « " 


(A-43) 


H  <<:  H  <<  T’* 


(A-46) 


Ln.- 


If  the  field  Is  written  as  in  Eq.  (A-43),  the  partial  polarization  may  be 
written  in  terms  of  in-phase  and  out-of -phase  components  C  and  S,  respec¬ 
tively,  as 

P(k,z,t)  ■  C(k,z,t)  cos(vt  -  nz  +  $)  +  S(k,z,t)  sin(vt  -  nz  +  $)  (A-47) 


-.v 
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o 


(A-54) 


where  cq  Is  Che  permittivity  of  free  space. 

The  equilibrium  values  of  the  density  of  states  and  p are 
determined  by  the  quasi-Fermi  levels.  The  Fermi  levels  are  determined  by 
temperature  and  the  voltages  maintained  across  Che  device.  Expressions  for 
these  equilibrium  densities  are 
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(A-55) 


m 


fc(1  -  fv> 


CO  2 

IT 


V1  -  fc> 


VO  2 

it 


where  fc  and  fy  are  the  Fermi  functions  shown  in  Fig.  A-3 


(A-56) 


f 

c 


1 


2  2 

1  +  exp[(Eg  -  Fc)/KT] 


(A-37) 


f 

v 


2  2 

1  +  exp[(i^  -  FJ/KTJ 


(A-58) 


We  take  the  quantity  KT  as  1/40  eV.  Eq  is  the  band-gap  energy  of  the 
intrinsic  region  and  is  1.6734  eV,  the  quasi-Fermi  level  for  the  conduction 
band  is  0.0172  eV.  and  the  quasi-Fermi  level  for  the  valence  band  is  1.8536 
eV.  The  mass  of  the  electron  state  in  the  conduction  band  is  mc  =* 

7.61454  x  10~^g.  an<j  Che  electron  mass  in  the  valence  band  niy  * 

4.9367  x  10~^8g.  As  long  as  the  voltage  and  temperature  are  maintained,  the 
equilibrium  values  for  the  population  will  remain  constant.  This  issue  will 
be  addressed  again  later  in  the  report. 

At  this  point  we  introduce  retarded  time  "x,"  transforming  to  a  set  of 
new  variables 


T  " 


(A-59) 


z  ■  z 


(A-60) 


so  that 
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- -  -T*  n  C  —  *  — 

3z  o  3 1  3  z 


3_  3_ 

3t  *  3t 


and  we  have 


3<f 

32 


1/2  n 


/ 


n  e 
o  o 


*[&]  -  -i/2  -/  ~r~ 

n  e 
o  o 


3S 

3t 


"  -YCVS  -  (u>  -  v  -  *)C  -  ^(Pc  -  Pv) 

If  ■  -YcvC  +  («  '  V  -  *)S 
3p 

FT  ■  'Wc  -  ‘’c’  •  Vc  +  i  T 
FT  *  -Tcol(,,v  •  pvo>  +  Vc  •  I  T 


(A-61) 


(A-62) 


(A-63) 


(A-64) 


(A-65) 


(A-66) 


(A-67) 


Generally,  Eqs.  (A-62)-(A-67)  would  be  the  set  of  equations  used  for  calcula¬ 
tions.  However,  because  of  the  size  of  y  and  y  . ,  further  modification  and 
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approximations  are  necessary.  Since  we  will  be  interested  in  pulses  whose 
widths  are  generally  greater  than  1  psec  and  Ycv  is  of  the  order  of  101-*,  we 
can  neglect  the  time  derivatives  of  S  and  C  in  Eqs.  (A-64)  and  (A-65). 

Further  $  can  be  dropped  from  these  equations.  This  immediately  implies  that 
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m 


p2 

s  -  - _ , _ 21 _ _  (p  -  p  ) 

"  +  (»  -  V)2]  c  v 


C  -  -  P2  ■?(“-  V)  ,r  .  , 

t[v  /  +  (.  -  ZF\  c  v 


(A-68 ) 


(A-69 ) 


One  might  assume  that  the  time  derivatives  In  Eqs.  (A-66)  and  (A-67 )  could 
also  be  set  to  zero  within  the  same  approximation.  This  would  Imply  that  the 
density  of  states  for  the  upper  and  lower  states  was  essentially  given 
by  pco  and  Pv£).  However,  we  assumed  Initially  that  PcQ  and  Pyo  were  deter¬ 
mined  by  the  voltage  across  the  intrinsic  region  and  that  this  voltage  remains 
constant.  In  reality,  as  the  population  levels  are  drawn  down  because  of 
stimulated  emission,  the  voltage  across  the  region  changes.  Inducing  a  current 
flow  in  the  device.  The  intrinsic  capacitance,  resistance,  and  mobility,  as 
well  as  external  circuit  parameters,  determine  how  the  Fermi  levels  are 
restored  to  the  equilibrium  value  of  Pcq  and  PyQ.  This  would  entail  the 
addition  of  extra  equations  to  take  these  effects  into  account.  As  an 
approximation,  we  will  replace  all  three  extra  equations  and  Eqs.  (A-66)  and 
(A-67)  by  simply  changing  Ycq1  to  yp,  where  yp  is  an  equivalent  pumping  rate 

to  the  system.  We  take  y  ■  5  x  109/sec.*  This  gives  the  following  set  for 

P 

equations: 


.  ,  dkycv(Pc  ~  Pv) 

«  -L  +  — * —  /  - - - — 

2  n  e  [y  +  (w  -  v)  ] 
o  o  1  cv  J 


3<fr  JW2S  r  dk(u  ~  v)(Pc  ~  Pv) 

w"5 K77J  [y  2+(„-.)2] 

o  o  1  cv  J 


3PC  ~  pv) 

7T-  '  -VPc  -  Pco>  '  Vc  *  2  +  (B  .  v)*j 


(A-70) 


(A— 7 1 ) 


(A-72) 


H.  Haus,  private  communication. 
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(A-73) 


where  we  have  replaced  the  symbol  p  by  y  and  assumed  that  this  dipole 
moment  y  is  independent  of  momentum  It.  The  limits  of  integration  for  Eqs. 
(A-70)-(A-7l)  and  the  value  of  y  will  be  discussed  later.  The  quantity  L 
represents  the  linear  loss  term  and  was  represented  by  k/no  earlier.  Two 
regimes  of  operation — steady  state  and  pulsed — are  of  interest.  These  two 
regimes  will  be  considered  in  the  following  sections. 


4.  STEADY-STATE  OPERATION 

For  steady-state  operation,  the  time  derivatives  of  p  are  set  to  zero, 

and  the  population  difference  N  *  p  -  p  is  immediately  given  by 

c  v 

- — 7 -  P  “  P 

Y  +  Y  CO  VO 


i  +  — l€l  2 

*  \v(ys  +  y 


where 


*U>)  - 


Ycv2  +  (w  '  v)2 


(A-74) 
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By  defining  a  dimensionless  Intensity 


*2v»<'r. +  V 


Eqs.  (A-70)  and  (A-71)  can  be  rewritten  as 
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The  factor  k2  is  shown  explicitly  in  the  integrals  of  Eqs.  (A-77)  and  (A-78). 
The  value  of  the  dipole  moment  y  may  be  determined  experimentally  in  an 
indirect  fashion.  It  has  been  determined  (Ref.  A-7)  that  the  quantity 
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From  the  values  of  reduced  mass  and  energy  of  transition  (1.51  eV)  used  in 
that  reference  and  an  index  of  refraction  of  3.62,  the  dipole  moment  is  u  * 
5.62  x  10~27  coulomb  -cm,  which  is  equivalent  to  a  separation  distance  of 
3.5  A.  Finally,  note  that,  if  the  limits  of  integration  are  allowed  to  go  to 
infinity,  then  the  integral  in  Eq.  (A-78)  is  unbounded.  An  expanded  dis¬ 
cussion  of  this  point  has  been  included  in  the  main  body  of  the  text. 

5.  PULSED  OPERATION 

Examination  of  Eqs.  (A-70)-(A-73)  leads  one  to  make  several  modifi¬ 
cations.  For  the  steady  state,  we  define  the  population  difference 
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and  also  the  sum 
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By  adding  Eqs.  (A-72)  and  (A-73),  we  obtain 
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where  we  have  replaced  t  by  t 
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Since  it  is  assumed  that  a  significant  time  will  have  elapsed  between  turning 
the  device  on  and  the  beginning  of  optical  amplification,  Q  will  have  reached 
its  equilibrium  value  Q0,  and  N  will  have  an  initial  value  equal  to 
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For  this  reason,  during  the  operating  time,  Eqs.  (A-72)  and  (A-73)  are 
identical  and  equal  to 
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where  Nc  ■  pcq  -  pv(j.  The  population  difference  then  has  a  formal  solution 
given  by 
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Between  pulses  when  the  field  is  zero  (since  we  have  assumed  that  the  polari¬ 
zation  follows  the  field  exactly),  the  population  density  will  have  a  time 
dependence  given  by 
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where  Nint  in  Eq.  (A-88)  or  Nint(tcP  in  Eq.  (A-89)  is  the  population  differ¬ 
ence  at  the  beginning  of  a  pulse.  It  could  also  be  the  difference  just  at  the 
end  of  a  pulse  where  the  population  again  grows  according  to  Eq.  (A-89). 

This  allows  one  to  determine  the  population  density  at  the  beginning  of 
each  pulse  without  having  to  perform  integrations  between  pulses.  If  the 
complex  field  S  is  defined  as 
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then  Eqs.  (A-70)  and  (A— 7 1)  are  rewritten  as 
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For  pulsed  operation,  Eq.  (A-91)  is  solved  for  the  real  and  Imaginary  parts 
given  the  results  of  Eq.  (A-88)  for  N(k,t,z).  The  numerical  procedures  used 
to  solve  Eq.  (A-91)  differed  considerably  from  the  technique  used  for  the 
steady-state  case.  A  hlgh-order  Runge  Kutta  was  used  to  perform  the  z  Inte¬ 
gration.  A  straightforward  trapezoidal  rule  was  used  to  perform  the  time 
Integrations  in  Eq.  (A-88).  The  momentum  Integration  was  performed  by 
assuming  a  polynomial  fit  to  the  third  order  in  the  k  dependence  of 
N(k,t,z).  This  third-order  polynomial  was  then  integrated  analytically.  The 
steps  In  the  process  Involved  are  that  at  each  z  step  a  set  of  values  for  the 
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electromagnetic  field  at  all  given  values  of  time  for  some  previous  z  value 
are  prescribed.  These  values  of  <f  are  then  used  to  determine  the  population 
difference  at  the  given  values  of  time  and  momentum.  The  integral  over 
momentum  is  performed  and  the  new  set  of  values  are  determined.  Numerical 
results  will  be  given  after  a  discussion  of  device  considerations. 

6.  DEVICE  CONSIDERATIONS 

In  this  report,  we  have  used  various  values  for  the  band-gap  energy, 
electron  mass  in  the  different  bands,  and  Fermi  levels.  In  Section  2  and 
Fig.  A-2,  a  specific  model  of  a  three-layer  P-n-N  device  was  given.  The 
choice  of  this  specific  model  was  based  on  parameters  for  aluminum  concen¬ 
tration  and  doping  levels  for  devices  grown  during  the  experimental  inves¬ 
tigation,  which  were  carried  out  concurrently  with  this  work.  The  steps  taken 
to  obtain  the  various  parameters  given  are  taken  from  Ref.  A-l,  Chapter  4. 

First,  the  band-gap  energies  for  the  direct,  L,  and  X  bands  were  deter¬ 
mined  for  both  the  wide-  and  narrow-gap  materials.  Next,  the  electron  and 
hole  masses  were  determined  for  the  conduction  and  valence  bands.  Both  band- 
gap  energy  and  mass  are  functions  of  aluminum  concentration.  The  jumps  in 
conduction  and  valence  band  energies  were  taken  as  0.85  and  0.15  of  the  dif¬ 
ference  in  direct  band-gap  energies.  The  Fermi  levels  were  then  determined 
from  expressions  relating  the  electron  and  hole  carrier  concentrations  to  the 
temperature ,  carrier  mass,  band-gap  energies,  and  Fermi  level.  Assuming  con¬ 
stant  temperature,  the  Fermi  levels  can  be  determined  by  graphical  solution. 

In  the  intrinsic  region,  a  similar  procedure  gives  the  location  of  the  Fermi 
level,  though  for  this  region  the  exponential  approximation  to  the  Fermi 
integrals  may  be  used.  The  built-in  fields  and  potentials  are  determined  from 
Poisson's  equation — one  for  each  region.  It  is  assumed  that  the  mobile  car¬ 
riers  can  be  neglected  in  the  space  charge  region.  By  making  use  of  the  con¬ 
tinuity  of  electric  flux  density  and  potential,  we  can  express  the  potential 
in  each  region  of  the  P-n-N  device  in  terms  of  the  built-in  potential,  the 
widths  of  the  depletion  region,  the  dielectric  constants,  and  the  charge 
densities.  It  is  assumed  that  the  entire  Intrinsic  region  Is  a  space-charge 
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region.  By  using  the  continuity  of  the  electric  flux  density  at  the  junction 
of  the  n-N  region,  we  can  relate  the  depletion  widths  in  the  P  and  N  regions. 

Finally,  using  the  equality  of  the  total  built-in  potential  and  the  sum 
of  the  built-in  potential  in  each  region,  we  can  determine  the  depletion  width 
in  the  N  region  by  solving  a  quadratic  equation.  The  total  built-in  potential 
Is  1.836  eV,  which  in  our  case  is  the  difference  between  the  wide-gap  direct 
band-gap  energy  and  the  sum  of  the  Fermi-level  energies  of  0.085  and  0.064  eV 
(Eq.  4.3-8,  Ref.  A-l).  It  turns  out  that  the  device  looks  almost  like  a 
capacitor  with  most  of  the  voltage  drop  across  the  Intrinsic  region  (1.798 
eV).  The  built-in  potential  on  the  P  side  is  0.0193  eV  and  on  the  N  side  is 
0.0191  eV,  whereas  the  corresponding  depletion  widths  are  4.986  and  4.961  pm, 
respectively.  With  these  numbers  and  the  equations  for  the  potentials.  Fig. 
A-2(A)  was  drawn.  If  the  device  is  biased  by  1.836  V,  then  the  band  structure 
looks  like  Fig.  A-2(B).  It  is  assumed  that  this  is  the  usual  operating  con¬ 
dition.  The  assumptions  leading  to  the  band  structure  in  Fig.  A-2(A)  do  not 
hold  exactly  for  the  completely  biased  case  in  Fig.  A-2(B),  because  the 
intrinsic  region  is  not  an  exactly  uniform  depletion  region. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting  exper¬ 
imental  and  theoretical  Investigations  necessary  for  the  evaluation  end  applica¬ 
tion  of  scientific  advances  to  new  military  space  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  In 
dealing  with  the  many  problems  encountered  In  the  nation’s  rapidly  developing 
apace  systems.  Expertise  In  the  latest  scientific  developments  is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  snd  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  high- temperature  thenaomechanlcs ,  gas  kinetics  and  radiation;  research 
in  environmental  chemistry  and  contamination;  cw  and  pulsed  chemical  laser 
development  Including  chemical  kinetics,  spectroscopy,  optical  reaonatora  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  In  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bioenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  CaAs  low-noise  and 

power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  imaging;  mllllmeter-wave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation,  performance-sensitive  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  Intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  msterlala  In 
apace  environment;  materials  performance  in  apace  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  in  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  end  Ionospheric  physics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
waves  in  the  magnetosphere;  solar  physics,  infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  Ionosphere,  end  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  In  space  on  apace  systems. 


